arXiv:1504.04501vl [math-ph] 17 Apr 2015 


DCPT-15/17 


Algebraic Bethe ansatz for Q-operators: 
The Heisenberg spin chain 


Rouven Frassek 


rouven.frassekOdurham.ac.uk 


Department of Mathematical Sciences, Durham University, 
South Road, Durham DH1 3LE, United Kingdom 


Abstract 

We diagonalize Q-operators for rational homogeneous sl(2)-invariant Heisen¬ 
berg spin chains using the algebraic Bethe ansatz. After deriving the fundamen¬ 
tal commutation relations relevant for this case from the Yang-Baxter equation 
we demonstrate that the Q-operators act diagonally on the Bethe vectors if the 
Bethe equations are satisfied. In this way we provide a direct proof that the 
eigenvalues of the Q-operators studied here are given by Baxter’s Q-functions. 



1 Introduction 


Integrable spin chains are prominent examples of integrable models. The solution of the 
Heisenberg XXX 1 spin chain goes back to Hans Bethe [J]. His ansatz of factorized plane 
wave scattering is today known as the coordinate Bethe ansatz and is probably the most 
intuitive way to solve the eigenvalue problem. However, the origin of the integrable struc¬ 
ture remains unclear. 

The situation is somehow different when studying integrable models using the quantum 
inverse scattering method which, in particular, provides a prescription of constructing 
commuting families of operators. Combined with the algebraic Bethe ansatz, that allows 
to determine the wave function for models with a suitable reference state, it yields a 
powerful tool to diagonalize integrable Hamiltonians along with the commuting family of 
operators, see e.g. [2j. 

Q-operators are distinguished members of the family of commuting operators. Their 
eigenvalues, so-called Baxter’s Q-functions, are in general given by polynomials with zeros 
at the Bethe roots [3]. They are of fundamental importance for the so-called fermionic 
basis, see e.g. [3], as well as in the ODE/IM correspondence [5]. Furthermore, the spectral 
problem of planar J\f = 4 super Yang-Mills theory in four dimensions, which is believed to 
be exactly solvable in the planar limit [B], resulted in a set of relations among Q-functions, 
see and references therein. These all-loop equations are the same as for the spin chain, 
a close cousin of the Heisenberg chain studied here, that appears in the spectral problem 
at one-loop level. Initially Q-operators were introduced by Baxter to solve the eight-vertex 
model without specifying a reference state |3j. However, the limit to the rational case is 
singular. 

Inspired by [Sj, Q-operators were built in the framework of the quantum inverse scatter¬ 
ing method for finite-dimensional rational spin chains in [hi ITT] . The construction naturally 
incorporates a magnetic field (j) that manifests itself as a diagonal twist in the transfer 
matrices, see also ra- The twist field lifts the usual degeneracy in the spectrum and 
allows in the case of sl(2) for the definition of two Q-operators Q-t which belong to the 
commuting family of operators. A main ingredient in the construction mu are certain 
factorization formulas that relate the Q-operators to transfer matrices and yield the quan¬ 
tum Wronskian on an operatorial level. Having derived the quantum Wronskian and in 
particular the QQ-relations on the operatorial level it was argued that the eigenvalues of 
the constructed Q-operators are indeed the Q-functions, see Section 3.5 in [9]. Further¬ 
more, it was shown in [13] how local charges can be obtained from the Q-operators. For a 
discussion of Q-operators for supersymmetric spin chains relevant in the spectral problem, 
see [T31 - HB] . 

In the following, we diagonalize the Q-operators of the Heisenberg XXX S spin chain 
introduced in pin] directly using the algebraic Bethe ansatz. In analogy to the algebraic 
Bethe ansatz for the corresponding transfer matrix we define the off-shell Bethe vectors 
as the m-fold excitations on the tensor product of highest or lowest weight states. In order 
to evaluate the action of the Q-operators on the off-shell Bethe vectors we derive certain 
fundamental commutation relations from the Yang-Baxter equation. We show that the 
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Q-operators act diagonally on the Bethe vectors 


Q±(y)\^m) = Q±(yMm), 

with the spectral parameter y and Baxter’s Q-functions 

m 

Q±(y ) = e ±iy *H (y-zf), 
j=i 

if the Bethe roots zf satisfy the Bethe equations 

4 + s _ + 2 a fr 4 - zf + 1 

4 - s M 4 - 4 - 1 

This pragmatic approach motivates the choice of Lax operators used to construct the 
Q-operators, avoids the detour of deriving the functional relations from the factorization 
formulas and provides a transparent derivation of the eigenvalues of the Q-operators along 
with the corresponding eigenvectors. 

This article is organized as follows: In Section [2 we provide a short review of the Q- 
operator construction as presented in PHD- The elementary solutions used to construct 
the Q-operators satisfy certain Yang-Baxter relations. This is discussed in Section |HJ In 
particular, the Yang-Baxter equation relevant to derive the fundamental commutation 
relations in Section [I] is derived. Section 0 contains the algebraic Bethe ansatz for one of 
the Q-operators of the Heisenberg XXX S spin chain. We show that the Q-operator Q + 
acts diagonally on the Bethe states if the Bethe equations are satisfied. Finally, we use 
a relation between the Q-operators to obtain the eigenvalues and eigenfunctions of Q . 
Further material is collected in the appendices. 


(1.1) 

( 1 . 2 ) 

(1.3) 


2 Q-operators for the Heisenberg spin chain 

In this section we construct the Q-operators for the Heisenberg XXX S spin chain follow¬ 
ing pm], The elementary Lax operators used to construct Q-operators of the twisted 
Heisenberg spin chain can conveniently be written as 

TZs,±(y) = ^ aS ^l(y) e~ aS± , ( 2 . 1 ) 

with the oscillator algebra and the 51(2) generators obeying the commutation relations 

[a, a] = 1 , [S+, S_] = 2 S 3 , [S 3 , S ± ] = ±S± . (2.2) 

The diagonal part Tl±(y) depends on the spectral parameter y and reads 

= (2 - 3) 
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Here we have chosen a suitable normalization for highest/lowest weight representations, 
i.e. for 

5±|k ;± ) = 0, *S , 3|w ± ) = |oj = i= ) , (2.4) 

we have = |o; =1= ). This normalization ensures that all entries of the Q-operators 

will be polynomials in the spectral parameter when considering finite-dimensional rep¬ 
resentations s = 1, |,... up to an overall factor, see below. The Lax operators (12.11) 

are also employed in the so-called DST-chain, see e.g. mi- However, here the role of the 
auxiliary and quantum space is interchanged. The auxiliary space of the monodromy of 
the Q-operators for rational spin chains is built from the product in the oscillator space 
and from the IV-fold tensor product in the sl(2) space 

M±(y) = u s ,±{y) ® T^a,±{y) ® ® K a ,±(y) • (2.5) 

'-v-' 

N 

Here we used <g) to denote the tensor product in the first space and multiplication in the 
second space of TZ S ,±- Following PEU, we define the Q-operators as the regulated trace 
of the monodromy (12.51) 


Q±(y)=e ± ^ZZ 1 tr 


e ^2i<f>aa 


M±(») . 


( 2 . 6 ) 


where the trace sums over all states in the Hilbert space and we introduced the normal¬ 
ization 


Z± = tr 


^f2i(paa 


1 

l — e =F2 i<j> ' 


(2.7) 


The parameter is the twist angle and manifests itself as a diagonal twist in the transfer 
matrices. It is known that such diagonal twisted boundary conditions break the sl(2) 
symmetry but the quantum number m is still well defined. As a consequence, not only 
highest weight states but also all decendents can be obtained from the Bethe ansatz. 


3 Yang-Baxter relations 

The fundamental relation underlying the quantum inverse scattering method and the 
algebraic Bethe ansatz for rational spin chains is the Yang-Baxter equation. It reads 

Ai, 2 (x - y)R h 3 (x)R 2 , 3 (y) = R 2 ,3(y)Ri,3( x ) R iA x ~ y ) > (3- 1 ) 

and is defined on the tensor product of three vector spaces V 3 ® V 2 <S> V 3 . Here the R- 
matrix R hJ acts trivially as the identity on 14 with k ^ i,j and non-trivially on Vi®Vj. 
For further references we refer the reader to the excellent introduction [2]. 

The Lax operators C relevant to construct the transfer matrix of the Heisenberg chain 
with the fundamental representation in the auxiliary space satisfy the RTT-relation which 
serves as a defining relation of the Yangian, see e.g. [TS]. Namely, it holds that 

Rn,n(a; - y)C>u,a{x)£nM = /f n , s (|/)£n, s (^)Rn,n(^ - y) . (3.2) 
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Here the fundamental R-matrix reads 




/ x + 1 0 0 0 \ 

Oxl 0 

0 1x0 

V 0 00 x+1/ 


(3.3) 


and acts trivially on the third space denoted by s in (Id.21) . The sl(2)-invariant Lax operator 
can be written as 


£n, s (y) 


V + 1 + S3 

S+ 


s - 

f+i-s,J’ 


(3.4) 


with the generators of 51(2) defined in (12.21) . Here, the shift in the spectral parameter has 
been introduced for convenience. Its role will become clear later on. Before presenting the 
Yang-Baxter equation that is relevant to derive the fundamental commutation relations 
necessary to diagonalize the Q-operators in (12.61) we recall the derivation of the Lax 
operators TZ S ,± (EHj) employed in the construction [SJHTj. Apart from the solutions that 
are classified by the Yangian, the RTT-relation allows for solutions where the dependence 
on the spectral parameter is not proportional to the diagonal, cf. HD) as well as p3] where 
these solutions were studied in the context of Drinfcld’s second realization. I 11 particular, 
it holds that 


Ro.nO - y)Ln,±{x)La t ±(y) = An,± (?/)£□,±(x)R n , n (x - y ), (3.5) 

with 

= ( L Aaa ) “ d L °.-M = ( Z ~T t ' < 3 ' 6 > 

These solutions where used for the Q-operator construction in £0]. Furthermore, they are 
relevant to obtain the Lax operators for Q-operators for arbitrary representations of sl(2) 
in the quantum space in HP, see also Ha- Another ingredient is the sl(2)-invariant Lax 
operator for the transfer matrix of rational spin chains (13. 4p . The relevant Yang-Baxter 
relation that defines the Lax operators 1Z S) ± in (12.11) reads 


CnA x ~ y)La,±(x)U a>± (y) = K a ,±(y)Ln,±(x)£n, s (x - v) ■ ( 3 - 7 ) 

We note that for spin \ we identify IZi ± {z + |) = L u ,±(z). However, this relation is not 
suitable for our purposes. 

To obtain the fundamental commutation relations that are needed to diagonalize Q- 
operators we define the Lax operator 


£s,m(£) 


* + S3 S. \ 
S+ z — S 3 J 


It satisfies the unitarity relation 

c 

C,p(x - y)C n ,,(y - x) = (y - x)(x - y - 1) + -, 


(3.8) 


(3.9) 
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where c is the Casimir element of sl(2) given by c = S+ S- + S- S + + 2 S 3 S 3 . We note that 
the Lax operator (13.811 satisfying (13.9jl only differs by a shift in the spectral parameter 
from the Lax operator defined in (13.411 . This is not always the case and for convenience 
we distinguish between the two solutions (13.411 and (13.811 . In addition to the Lax operators 
L n ,± in (I3.6|l we define 


L+,n( z ) 


z T 1 T ua 
a 



and L- t n(z) 


that obey the unitarity conditions 


-1 a \ 

—a z + 1 + aa I ’ 


(3.10) 


L+,u(x-y)L a.+(y - x ) = x -V. L-,n(x-y)Lu,-(y- x )-x-y. (3.11) 

In contrast to the Lax operators C these solutions significantly differ from (13.611 . The Yang- 
Baxter equation relevant to derive the sought-after fundamental commutation relations 
can then be obtained from (13.7p after multiplying (13.811 and (13.1011 . We find 


'R' S ,±{y)C s ,n{x)L±,n{x - y) = L± ta (x - y)£ s p(x)n Sj ±(y) , (3.12) 

where we used the unitarity relations (13.911 and (13.1111 and for convenience renamed the 
spectral parameters. 


4 Fundamental commutation relations 


Following the algebraic Bethe ansatz we define the monodromy of the fundamental trans¬ 
fer matrix 0 as the matrix product of the Lax operators introduced in ( 13.811 in the two 
dimensional auxiliary space and the iV-fold tensor product in the quantum spaced la¬ 
belled by the representation label s. Using the same conventions as in the definition of 
the monodromy of the Q-operators ( 12.511 it can be written as 

M{y) = £ s ,n(y ) ® £ s ,a{y) ® • ■ ■ <8> £ s ,n(y) ■ (4-1) 


By construction Ai is a 2 x 2 matrix in the auxiliary space with operatorial entries built 
from the sl(2) generators at the different sites 

( My) B (y) \ 

\c(y) My) ) 

The corresponding transfer matrix is constructed as the trace over the two-dimensional 
auxiliary space containing a suitable boundary twist that acts trivially in the quantum 
space 

T n (y) = tr e~ 2i ^M(y) = e~^A(y) + e + ^D (y) . (4.3) 

Tn contrast to the convention often found in the literature where the fundamental transfer matrix 
refers to the transfer matrix of equal representation in quantum and auxiliary space. 

2 Note that the monodromy in (14.11) is related to one built from the Lax operator in (13.41) by a shift in 
the spectral parameter, cf. (1TSD . 



(4.2) 
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As a consequence of the Yang-Baxter relation (13.21) and the unitarity relation (13.91) the 
monodromy JV[ in (14.11) satisfies the RTT-relation 

R-n.nO - 2/)(I <8) M(x))(M(y) 0 1) = ( M(y ) <g> I)(10 M(x))R 0 ,n(x - y ), (4.4) 

where I denotes the 2x2 identity matrix. To obtain the fundamental commutation re¬ 
lations we insert the explicit form of the R-matrix Rn,n as given in (l3.3l) into (l4.4l) . see 
e.g. [2U2D]- In particular, we find 


B(x)B(y) = B(y)B(x), (4.5) 

A(x)B(y) = f(y, x ) B(y)A(x) - g(y, x ) B(x)A(y), (4.6) 

D(z)B( 2 /) = f(x, y) B(j/)D(x) - g(x, y) B(x)D(y), (4.7) 


C(a;)B(y) = B(y)C(x) +g(x,y) [A(y)D(x) - A(x)D(y)] . 


Here we used the notation 


f(x,y) 


—- and g(x,y) 
x-y 


1 

x-y ’ 


(4.8) 


(4.9) 


with f(x,y) = g(x,y ) + 1. 

To diagonalize the Q-operators in (12.61) we derive the commutation relations of the Q- 
operator’s monodromy M± with the operators A, B, C and D in (14.21) . As a consequence 
of the Yang-Baxter relation (13.121) . we hnd that the monodromy of the Q-operators (12.51) 
and the monodromy of the fundamental transfer matrix (14.11) satisfy 


M±(y) M(x) L±'U(x -y) = L±n(x - y)M(x) M±(y ). (4.10) 


In the following we focus on the Q-operator Q + built from the monodromy M + . In this 
case the additional fundamental commutation relations necessary to diagonalize the Q- 
operator Q + arise from (14.101) after inserting the explicit form of n as given in (13.101) . 
We obtain 


m 4 

-(y)l(x- 

■ y + 1 + aa)A(x) 

+ aB(x)] = [A(x)(a; 

— y T 1 T aa) — 

■C(x 

) a] M + 

(»), 








(4.11) 

M + 

(y)B(x) 

= B (x)(y — x — 1 

— aa) M + (y) + D(x) 

a M + (y) - M + 

(y) a 

A(x), 

(4.12) 

C(x 

)M +(y) 

= M+(y)(y-x- 

- 1 — aa)C(x) + A(x) 

aM + (y) - M + 

(y) a 

D(x), 

(4.13) 



[D(x), M + (y)] = B(x) aM + (y) + 

M+(y)aC(x). 



(4.14) 
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For our purposes it is enough to consider (j4.12ji and (14.141) and combine them into the 
relation 

M + (r/)B(a;) = (y - x)B(x) M+(y) + X(y, x ), (4.15) 

with 

X(y,x) = aM + (i/)D(i) + aM + (y) aC(x) — M + (y) aA(x). (4.16) 

The relations relevant to diagonalize Q involving M_ and L_p can be found in Ap¬ 
pendix [A] Similar relations appeared previously for the product of two Q-operators in ra- 

5 From Q-operators to Q-functions 

In this section we diagonalize the Q-operator Q + following the algebraic Bethe ansatz. 
In addition to the fundamental commutation relations (I4.5I) - (I4.8|) we will repeatedly use 
the commutation relation among the monodromy M + and the B-operators (14.151) that we 
obtained from the Yang-Baxter equation (14.101) . 

Following P], we introduce the reference state 

|D+) = |w + ) <g>... <8) |w + ), (5.1) 

N 

with the highest weight state at each site as defined in (12.41) . The action of the monodromy 
(14.11) on the reference state |f2 + ) yields in particular 

A(y)V + ) = <*(y)|W), D(t/)|n + ) = «(i/)|n + >, c(t/)|n+) = o, (5.2) 

with 

<x(y) = {y + s) N , S(y) = {y-s) N . (5.3) 

Due to our choice of normalization of the operators 7l± in (12 .3jl the action of the mon¬ 
odromy (12.51) of the Q-operator Q + on the reference state |f2 + ) is independent of the 
spectral parameter y and reads 

M + (j/)|n + )=e iS “|SJ + ), (5.4) 

where S l ot = J2iLi 5_ with S_ denoting the generator S_ acting on site i. The off-shell 
Bethe sates in the context of the algebraic Bethe ansatz are then defined as 

IV£> = B(z+)B(4) • • -B(4)|D + ). (5.5) 

The states |are eigenstates of the corresponding transfer matrix if the parameters 
Zi satisfy the Bethe equations [2]. Thus we expect that the Q-operators Q± also act 
diagonally as they belong to the family of commuting operators and obey the Baxter 
equation, see [S]. 
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Instead of studying the action of the Q-operator on the off-shell Bethe states (15.51) 
it is convenient to act with the monodromy and take the trace afterwards. Using the 
fundamental commutation relations in (14. 15j) we obtain 

m 

M +(y)l^m) = 11(2/ - U + ) b (a + ) m +(?/)I^ + ) + \y,i z+ }), (5-6) 

i= 1 


where 


m k —1 


Is/. {^ + }> = E n(»- z . + ) B (4)'"B(y_ I )X(y,y)B( 2 + +1 )...B(4)|SJ + >. (5,7) 


fc=l i =1 


Here, the operator X is defined in (14.l(ip and {z + } denotes the set of parameters zf ,..., zf. 
The role of the terms on the right-hand-side of (15.61) becomes clear after we have taken 
the trace over the auxiliary space and using the action of the monodromy M + on the 
reference state (15.41) . From the definition of the Q-operators (12.61) we find that 

m 

Q+(2/)l^m> = e + ^n (y-ztMm) +1 y,i z+ })<h ( 5 - 8 ) 

i= 1 


with 


\y,{z + })+ = e +i * l, Z?tT 


0 —2i(f)aa 


\v, U + }> 


(5.9) 


The term proportional to | iff) yields the Q-function Q + if the parameters zf are Bethe 
roots, see (11.21) . Thus, we expect that the unwanted terms \y, {z + })<t> vanish if the param¬ 
eters zf satisfy the Bethe equations (11.31) . 

From the commutation relations among the B-operators (14.5j) it is obvious that the 
left-hand-side of (15. 6 p is symmetric in the parameters z t as well as the first term on the 
right-hand-side. As a consequence also \y,{z + }) and \y, {^ + })</, given in (15.7p and (15.91) 
must be invariant under permutations of zf. This property will be important later. 

The evaluation of the unwanted terms (15.91) is more involved than in the algebraic 
Bethe ansatz for the corresponding transfer matrix. The reason is that while the operators 
A and D conserve the number of B-operators, the operators M + and C can annihilate 
one B-operator, see fl4.6j) - (l4.8p and (14.151) . After we have commuted all operators A, D, C 
and M + through the B-operators and evaluating them on the reference state as discussed 
(15.21) and (15.41) the unwanted terms in (I5.9p can be written as 


fk(y,4 *5™)) 


\y,U + })* = e +iy * E E hSE V'i — A y (t+1 >) • • ■ B(y (ro) ) St, |SV). (5.10) 

o-g<S m . k= 1 1 7 


Here we moved the oscillators multiplied to the left of the monodromy M + to the right 
by commuting them through the regulator e -2?9iaa and subsequently using the cyclicity of 
the trace. Finally, we can take the trace of M + on the reference state 


(g2 itf> _ 


sL l«+>- 


tre -x*Sa M + (t/)a t |Sl + ) 


(5.11) 



























Assuming that the states corresponding to different families {y + } are linearly indepen¬ 
dent, see also the discussion in m, we demand that all coefficients in (15.101) vanish 

F*iv, <(!),• ••■<„,) = 0 , ( 6 . 12 ) 

for all permutations cr e S m . In Section Hull we evaluate J-f for which the condition (15.121) 
yields the Bethe equations (1 1.3 1) and in Section I7T21 we show that all other Tf with k > 1 
vanish if the Bethe equations are satisfied. 

There exist exceptional solutions to the Bethe equations which in particular include 
coinciding Bethe roots or so called exact string solutions which may yield null vectors or 
singular solutions, see e.g. I'2T. 22 and references therein. To some extend the situation is 
more feasable in the case of twisted boundary conditions. However, in the following we 
focus on the generic case where zf zf and zf — zf ^ ±1. 

5.1 Bethe equations: The first term 

Having in mind that the unwanted terms in (15.71) are symmetric in the parameters zf 
we calculate the first term Tf employing an argument similar to the one used to obtain 
the unwanted terms for the corresponding spin chain transfer matrix, see MU- Let us 
consider an ordering of the parameters zf as in (15. 7|) and focus on the first term in the 
sum where k — 1. In particular, all unwanted terms (15.101) that are independent of B(z^) 
must arise from it after taking the trace! We denote the equality of the part independent 
of B(zi“) as ~ such that 

\y,{ z+ }) -x(y,zf)\i) with l J ) = II 5 (4)l fi )> ( 5 - 13 ) 

fee/ 

as discussed above. Here the set / is given by / = (2, 3,..., m}. 

Let us evaluate the action of the operators A, D and C on |J) in X(y, zf)\I), cf. (14.161) . 
This can be done by repeatedly applying the fundamental commutation relations (14.6!) . 
(14.71) and (14.15j) and evaluate the relevant operators on the reference state (15.21) . The 
action of the operators A, D and C on |J) is well-known, see e.g. [2D], and can be found 
in (IB.2D (1B.4I) . We obtain 

X(y, y )|/> = (v'(zt) a M+(y) - A'tf) M + (y) a) |/> 

+ £ c /(*i>M + (y)a|/\{j}) 

jei 

+ E (V'M) aM + - A'(zt) M+ a) |(J \ {j}) U {1}) (5.14) 

3 ei 

+ E d*(y)aM + a|(/\{j,U)U{l}). 
j,kei 
j<k 

Here the caligraphic coefficients can be found in (IB .51) — (IB . 8D . As we can see, only the 
first term can contribute to J-f as all others already contain less than m— 1 B-operators. 
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When substituting A 1 and V 1 as given in (IB. 51) and (IB.61) into (15. 14j) we obtain 


(V'(zt) aM + (») - A'(4) M + (y) a) |J) = ]}(!/- 4W4) 

i =2 

M M 

X 5 (-i + ) I[f(4,4) a ^-(y) - a(zf) JJ f(zt,zf)M_(y)a 


1=2 


i =2 


|fi) + ■ ■ • ■ 
(5.15) 


Here the dots denote terms that contribute to Ajf with k > 1 and as such contain less 
B-operators. Finally, we evaluate the trace over the auxiliary oscillator space, cf. (15.8p and 
(12.61) . on the reference state (15.111) and obtain 


(v, z t, ••■,4) = e + ^ n(j/ - 4 ) \e +t *K4 ) II f(4,4)~e 


i =2 




44) Ilf (4,4 


i =2 


i=2 


(5.16) 

We demand that the coefficients vanish for all permutations of the parameters z+, 
cf. (15.121) . This yields the Bethe equations 


Q44 ~ l )44) + Q44 + for k = l,...,m, (5.17) 


compare (ll.3l) and (l5.3l) . Here we rewrote the product of functions / in (l4.9p for zf ^ z£ 
in terms of Baxter’s Q-function as given in (II.2p . 


5.2 Ravelled Bethe equations: Recursion 

Naturally, we expect that also the coefficients Jqf with k > 1 vanish if the Bethe equations 
are satisfied. This is indeed the case. However, the coefficients are quite involved and we 
will not present them here. In order to show that they vanish we again consider the first 
term in the sum on the right-hand-side of (15.151) and isolate the part that does not depend 
on B(^). As B^^) appears on the right of the monodromy M + in the last two terms of 
(15.141) we have to keep these terms and evaluate them further using (15.61) . We fold 

X(y,zt)\I) ^ (D J (4)aM + ( ! /)-/(4)M + ( ! ,)a) |/> 

+ 51 4 ( 4 ) a M +(y) a \4 0'}) 

jei 

+ Y, (Pj(4) 4) - 4(4) x (y, 4 ) a ) I AO'}) ( 5 - 18 ) 

iei 

+ 55 4A4) aX i.y,4) a \i\44 }), 

j,k£l 

j<k 


where we dropped the terms with B(^) on the left of the monodromy M + . From (15.181) 
and (I5.13P it becomes clear that the part which does not depend on B(z+) of \y, {T + }) 
can be written as 


\y,{z + }) ~X(y,zt)\I) ~ 55 55 

rrGSm -1 k =1 


1 

|</fc|! |7 \ Jfc|! 


G Uj k )(y^4) W(!\Jk)) , 


(5.19) 
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where Jk = {2,..., k}. The coefficients Gj k are symmetric in zf for i G I/ Jk and Zj~ for 
j G Jk and depend on the monodromy M + and the oscillators a but not on the operators 
A, B, C and D. They are determined in the following. Substituting (15.191) into (15.181) we 
obtain the recursion relation 


G j k (y, 4) = J2 {4 ( 4 ) a g 4\{J}(?/. 4 ) - 4(4) G (y, 4 ) 

1=2 

+ E c i,j (4 ) a (y, zt) a, 

2<i<j<k 

for all permutations of the elements 2,..., m with the initial conditions 

G v(y, 4 ) = '^ 1 ( 4 ) aM +(y) -^(^M+^a, 


B/'-v+'i nAlil 


(5.20) 


(5.21) 


Gfofa, zt) = v‘ (z+)®AW( z +) a 2 M+ (h) + A'(z+)A™(z+) M+(y) a 2 . (5.22) 

Note that in (I5.22|) we expressed Cj in terms of A and V using (IB.71) . The operators G j k 
are obtained recursively in Appendix O If we substitute A and V as given in (IB. 51) and 
(IB.61) they read 


Gj(y,z+,zj) = \J\lU44At) 

ieJ 


a i j i +i M + (y) n 5(z+) n f(z+, z +) 

leJu{i} iei\J 

-M + (?/)a |j|+i n 44 ) n f(4,4) 


jeJ u{i} 


iei\J 


(5.23) 


In principle, one can evaluate the coefficients J-)f by subsequently applying (15.191) and 
taking the trace as in (l5.9l) . However, this is not necessary! Taking the trace of (15.191) and 
moving all oscillators to the right of the monodromy M + yields 

m pi(y+k)4> y\k n ( 7 +, s -+) 

|y ’ {z+}) * “X,£— (524) 

x tr [ e - 2 ^. M+ to)a‘] B« (t+1) ) ■ • • B( z + m) )|!J+). 

with the coefficients 

km km 

Qt(y,zf,4) = e+ik * II 5 (4) n f(4^4)-^ ik<f> U44) A 44,4)- 

j =1 l=k +1 j =1 l=k +1 

(5.25) 

We note that up to a factor Qf coincides with and thus vanishes. To show that all 
with k > 1 vanish we note that the Bethe equations (15.171) imply 


,+<m ■ 


II S (4 ) 144 - 4 + !) - e n 44 ) IM - 4 - !) = 0 1 (5-26) 


16/ i e i 


16 / 


lei 
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where we introduced the set / C {1,..., m} of cardinality |/| and its complement I. From 
(15.261) it follows that 

GkiVi zi, 2 ), • • •, z a{m) ) = 0 . (5.27) 

As a consequence all unwanted terms that are independent of B(z^) vanish, cf. (15.241) . 

We discussed earlier that the unwanted terms and in particular (15.71) are invariant 
under permutations of the Bethe roots. Following the same strategy presented here for 
the case not containing B^f) we can equally well evaluate the terms that do not contain 
B(z 2 ~), B(^ 3 ~), etc. The corresponding coefficients Qf. are permutations of (15.251) involving 
zf and thus vanish due to (15.261) . Finally, we conclude that in general the unwanted terms 
(15.101) vanish if the Bethe equations (15.171) are satisfied. 

6 The other Q 

To diagonalize the Q-operator Q we can in principle repeat the procedure explained 
in Section 0 using the commutation relations in Appendix [A] The natural choice of the 
reference state is such that it is annihilated by the B-operator and excitations are created 
by acting with C-operators, see below. However, for finite-dimensional representations 
we can make use of the spin-flip symmetry of the Q-operators. First, we note that Lax 
operators H s ,± can be related to each other by the transformation 

Kn a , ± {y)K~ x = n s ^(y), (6.1) 

defined as 

K S± K - 1 = , K S 3 K~ x = -S 3 . (6.2) 

As a consequence of (16.11) and the definition of the Q-operators in (|2.6j) we see that Q± 
are related through the transformation K acting on each site of the spin chain and a sign 
change of the twist held f 

KQ ± (?/)K _1 = Q T (y)\ ( p^_ ( j ) with K = K ® I< ®® I< . (6.3) 

N 

As shown in Section ED the Q-operator Q + acts diagonally on the Bethe states \ if the 
Bethe equations (15.171) are satisfied. As a consequence of (16.31) we obtain the action of the 
Q-operator on the Bethe states 

W m ) = C(zf)...C(z-)\n~), (6.4) 

with |H _ ) = \oj~ <g)... <g) |a; - ) = K |H + ). The Q-operator Q_ acts diagonally on the Bethe 
states (E3D as 

m 

Q-(S/)I 'I’m) = e ~ W(P II (2/ - Z k)\^m) > (6-5) 

k =1 

if the Bethe roots zf satisfy the Bethe equations 

Q-( z k ~ !)«(^) + Q-( z k + = 0, (6.6) 
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with a and 6 given in (15.31) . cf. (II.ID — (jl.31) . Here we used that the transformation K acts 
on the quantum space of the monodromy A4 in (14.11) as 


( My) B(y) \ 1 _ / D(y) C(y) 
1, C(y) D(y) J * { B(y) A(y) 


(6.7) 


7 Conclusion 

We diagonalized the Q-operators for the XXX S Heisenberg spin chain using the algebraic 
Bethe ansatz. In doing so we demonstrate that the eigenvalues of the Q-operators are 
given by Baxter’s Q-functions. In analogy to the algebraic Bethe ansatz for the corre¬ 
sponding transfer matrix, the results obtained rely on certain fundamental commutation 
relations that were derived from the Yang-Baxter equation and the cancellation of the 
so-called unwanted terms. Our findings consolidate the construction of the Q-operators 
as carried out in pm. Here the functional form of the Q-operators has been deduced 
rather indirectly from their functional relations. 

We have seen that the commutation relations derived here are convenient to diagonalize 
the Q-operators Q-t on the states \M) respectively as the eigenvalues are manifestly 
polynomials. However, it is not apparent that the eigenvalue will be a polynomials when 
acting with Q± on the states \^) respectively. This yields a curious relation among the 
Bethe roots as presented in Appendix [E], see also [121123]. It may be worth to study these 
relations in the zero twist limit and in particular at half-filling. 

Furthermore, it would be interesting to generalize our procedure to the case of sl(n)- 
invariant spin chains and their supersymmetric relatives mnm where the nested Bethe 
ansatz has to be applied. 

Finally, we hope that our proof carries over to the g-deformed case where a similar 
procedure has been studied in Bi- 
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Appendices 

A Fundamental commutation relations for Q_ 


The fundamental commutation relations arising from the Yang-Baxter equation (14.1011 for 
the Q-operator Q read 

M_(?/)[(x — y + 1 + aa)D(x) + a C(x)] = [D(x)(x — y + 1 + aa) — B(x) a]M_(y), 

(A.l) 

M_(y)C(x) = C (x)(y — x — 1 — aa)M_(y) + A(x) aM_(j/) — M_(t/) aD(x), (A.2) 
B(Y)ML(?/) = M -(y){y — x — 1 — aa)B(x) + D(Y) aM„(r/) — M_(y) aA(i), (A.3) 


[A(x), M_ (j/)] = C(x) a (y) + M_ (y) a B(x). (A.4) 

Combining (1A.2|) and (IA.4D yields 

M_(r/)C(a;) = (y — x)C(x)M_(y) + aM_(y)A(x) + aM_(y)aB(x) — M_(y)aD(x). (A.5) 


B Commutation relations 


In this appendix we calculate the action of the operators A, D and C on a given off-shell 
state 

|/> = IlB(**)|fi>. (B.l) 


kei 


We are mostly interested in the case where the set I is given by I = {2, 3,..., m}. One 
finds 


AMI/} = + £ AY) (I \ {/}) U {!}) , 

jei 


(B.2) 


D(^i)|/) = V , (z 1 )\T) + £©'M (AM) U {!}), (B.3) 

i€i 


CM|/> = EC'MIAU}>+ EdMOKAfeUMU}). (B.4) 

j&I j,k£l 

j<k 

The coefficients A, V and C appearing above read 

A^zx) = a(z 1 )\{f(z k ,z l ) , A I j (zi) = a(z j )g(zi,z j ) fj f(z k ,Zj ), (B.5) 

kei k£l\{j} 
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^Ei) = 5(zi) n f( z i> z k), 'D I j(zi) = 8(zj)g(z j ,z 1 ) n f( z i, z k), (B.6) 

k£l kel\{j} 

and 

C'(zi) = Aj(z i) I>' v (zi) + ®/(zi) A^(z i), (B.7) 

C'^i) = -4'( Zl ) ©T’fe) + X>'(*i) A^M . (B-8) 

The formulas above can be deduced from the commutation relations in (14.611 and (14.711 . 
see also [2D] for a more detailed explanation. 


C Proof of the recursion relation 

Let us first note that (15. 23j) is compatible with the initial conditions (I5.2ip and (15.2211 . In 
order to show that G j (l5.23ll solves the recursion relation (l5.2(Jll we write (l5.23ll as 



Gjiv-i 4 

) = 

)-a 7 j ( 2/ ,4). 

(C.l) 

we defined 






= \J\\W44,4) 

n 

1 J] f(z+,z+)4 J \ +1 M + (y), 

(C.2) 


ieJ 

JGJU{1} 

iei\J 

a j(». 4) 

= M ! II 9(z?,zt) 

n 44 '. 

) f] f(zf, Z +)M + (y)a.W +1 . 

(C.3) 


ieJ 

j eJu{i} 

iei\J 


and 


Our proof is based on two identities 

Gj k (y, zf ) = D (y, zf ) - Aj fc (y, zf ) 

= E iPj (■ 4 ) a D E\E (V, 4 ) + 4 (z +) A j^ j} (y, zf ) a 


3 =2 


(C.4) 


and 


E c lj ( z i ) a g 2\{5} (V’ z i) a = E ( v )(4 ) a Klxi} (j/» 4 ) 

2 <i<j<k j= 2 

+ E(E) D M$} E) a 


(C.5) 


Substituting (1C. 1 jl for GEE bi the recursion relation (15.2011 we see that if (1C.411 and 
(1C.5[) hold also (15.2011 is satisfied. 

To show (IC.4p and (1C.5ft we note that 


D j{y, 4) = + E v j (4 ) a D Eli E) > 

j£J 


(C.6) 
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A j (y, 4 ) = - E A j ( z t ) A Ah] (y, 4)*, 


(C.7) 


j&j 


with the coefficients A and T> given in (IB. 51) and (IB.61) . These relations rely on the curious 
identity 

P P 

£.ilf(x i ,x k )=p, (C.8) 

k =1 i^k 

which generalizes to the trigonometric case and can be shown by induction in p G N and 
using the partial fraction decomposition 


j X Xi 1 ^ ^ 1 j—j- Xk X{ 1 


i— 1 




X Xk Xk Xi 


(C.9) 


It immediately follows that (1C.4|) holds. Furthermore, we note that the term on the left 
hand side of (1C.51) can be written as 


E C lj(4 ) aG M{3 } (?/’E)a = E p .7 V (-i) aD M{3}(^ffi + ) 

2 <i<j<k i =2 j =2 

l¥=i 

- E v l (4 E A ^\4 A S\{i]}(s/> E) a, 

i=2 J=2 


(C.10) 


see Appendix [B] Again using (1C.61) and (1C.71) we arrive at (1C.51) which concludes the 
proof. 


D Hamiltonian 


For completeness we give the Hamiltonian for the corresponding XXXi Heisenberg spin 
chain. It can be obtained from the transfer matrix in (14.dp by taking the logarithmic 
derivative at the shift-point, see e.g. [23 • It reads 


N 




i+l) > 


(D.l) 


i= 1 


with the permutation operator P iji+ i = (I + S^ l+1 ^)/2 and the twisted boundary 

conditions 

_ . . . ~(N\ 

(D.2) 


-p _ —2 i(j}Si N - -p +2i</>si N ^ 
r n,N+ 1 — e 3 ryd 3 


Here the generators S W = a t site % are given in terms of the Pauli matrices a. 


16 





















E Relations between both sides of the equator 

Solutions beyond the equator were studied in [23]. Here we present a relation between the 
roots zf and zj as also obtained in [T2]J for s = 1/2. The derivation relies on the property 
that Q ± have common eigenvectors and the fundamental commutation relations (14.111) 
and (14.131) . For s = 1/2 we find 


OO 


E e~ 2i * q (y - q 

9=0 


1 \N TT _ (V Z i ) _ 

2 iAiv-^-q^y-zr-l-q) 


N—m 


II (y~4)' 


(E.l) 


We can set y = z£ and obtain the relations 




(4 - z { ) 


q =o 


=i (4 -Zi — q) (4 — z i - 1 - v) 


= 0 for k = 1,..., N — m. 

(E.2) 

This formula has been checked for a few examples. Similar relations hold for finite¬ 
dimensional representations of higher spin. One finds 


00 m (v — 

E e- 2 ^ M «- (y, q) II 7- _ % l) - , -t 

4o t=i (y- z i-q)(y- z i- 1 ~q) 


2 sN—m 


= II 


(E.3) 


2=1 


with the action of the diagonal elements of the monodromy on the reference state |Q ) 
given by 

(y, 5, q) = ^ 2 Fi(-2s, -q,l-s-y, 1)^ , (E.4) 

which can be expressed in terms of T-functions. 
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